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Abstract 


Velázquez, J. A., Troin, M., € Caya, D. (September-October, 
2015). Hydrological Modeling of the Tampaon River in the 
Context of Climate Change. Water Technology and Sciences (in 
Spanish), 6(5), 17-30. 


This work compares the hydrological modeling of the 
Tampaon River Basin (in east-central Mexico) with two 
hydrological models (SWAT and GR4J) and then evaluates 
the impact of climate change on the water balance of the 
basin. The calibration and validation of the models (over 
14-year periods) show that both performed satisfactorily 
when simulating daily flows. The results indicate that 
SWAT more precisely reproduces observed mean monthly 
streamflow while GR4J] overestimates it during the dry 
season and underestimates it during the rainy season. The 
analysis of the impact of climate change was performed 
by using climate ensemble simulations derived from the 
Canadian Global Climate Model (CGCM3) downscaled by 
the Canadian Regional Climate Model (CRCM). The climate 
simulations (after bias correction) were used as input data 
for both hydrological models for two periods: a reference 
period (1971 — 2000) and a future period (2041-2070). The 
results indicate a significant decrease in mean monthly 
streamflow in the Tampaon River Basin for the future period 
(-36 to -55%), as well as a decrease in maximum monthly 
streamflow (-34 to -60%) and minimum monthly streamflow 
(-36 to -49%). The results from this study provide an overall 
perspective of the potential impact of climate change on the 
hydrological response of the Tampaon River Basin. 


Keywords: Canadian Regional Climate Model, SWAT, GR4J, 
climate change impacts. 


Resumen 


Velázquez, J. A., Troin, M., € Caya, D. (septiembre-octubre, 
2015). Modelación hidrológica del río Tampaón en el contexto del 
cambio climático. Tecnología y Ciencias del Agua, 6(5), 17-30. 


Este trabajo presenta, en un primer paso, un estudio de 
comparación de la modelación hidrológica de la cuenca del río 
Tampaón (localizada en el centro-este de México) con los modelos 
hidrológicos SWAT y GR4J, y en un segundo paso, la evaluación 
del impacto del cambio climático en el balance hídrico de la cuenca. 
La calibración y validación de los modelos (en periodos de 14 años) 
mostró un desempeño satisfactorio de ambos en la simulación de 
caudales diarios. Los resultados mostraron que SWAT reproduce de 
manera más precisa el caudal medio mensual observado, mientras 
que GR4] lo sobrestima en la temporada seca y lo subestima en la 
húmeda. El análisis del impacto del cambio climático se realizó a 
partir de simulaciones provenientes del Modelo Climático Mundial 
Canadiense (CGCM3) regionalizado con el Modelo Regional 
Canadiense (CRCM). Las simulaciones climáticas (una vez 
corregido el sesgo) se usaron como datos de entrada a los modelos 
hidrológicos para dos periodos: uno referente (1971-2000) y otro 
futuro (2041-2070). Los resultados estiman una disminución 
importante del caudal medio mensual en la cuenca del río Tampaón 
para el periodo futuro (entre -36% y -55%). Además, se valora 
una disminución del caudal máximo mensual (entre -34% y 
-60%) y mínimo mensual (entre -36% y -49%). Los resultados de 
este estudio aportan una visión general del potencial impacto del 
cambio climático en la respuesta hidrológica de la cuenca del río 
Tampaón. 


Palabras clave: modelo climático regional canadiense, SWAT, 
GR4J, impacto del cambio climático. 
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Introduction 


Assessing the impacts of climate change on 
water resources is one of the major challenges 
of the 21* century. The expected changes in tem- 
perature and precipitation will have an impact 
on the hydrological cycle and, therefore, on 
water availability for human uses and economic 
activities. Mexico is particularly vulnerable to 
climate change, since many Mexican regions are 
sensitive to extreme climate events (Arreguín- 
Cortés $: López-Pérez, 2013). In addition, water 
demand for agricultural, urban, and industrial 
sectors in Mexico has grown much more quickly 
because of the rapid expansion of these sectors 
in recent decades (Mundo éz Martínez-Austria, 
1993), and climate change will make responding 
to this increasing demand for water more diffi- 
cult (Magaña éz Conde, 2000). 

Climate change projections for Central 
America (Mexico included) suggest increases 
in temperature, and increases or decreases in 
precipitation that will have an impact in fu- 
ture streamflow and water availability in the 
region (IPCC, 2014). The projected change of 
temperature depends on the considered emis- 
sion scenario and geographical location. For 
instance, Martínez-Austria (2007) presented 
an analysis of projected temperature in Mexico 
based on GCM outputs: in scenario A2 (high), 
the projected temperature increase is expected 
to be between 4*C and 6*C, while in scenario 
B2 (medium) it will be between 2*C and 4*C. 
Martínez-Austria (2007) also discussed the re- 
sults obtained by Morales, Magaña, Barrera and 
Pérez (2001); this latter study identified regions 
with the highest and the lowest projected incre- 
ment of summer temperature (i.e., Northwestern 
Mexico and the Chiapas-Tabasco Zone, and the 
Yucatan Peninsula, respectively). 

Despite Mexico's vulnerability to the po- 
tential impacts of climate change, few studies 
have evaluated the impact of climate change 
on the country's water resources. For instance, 
Mendoza, Villanueva and Adem (1997) divided 
the Mexican territory into twelve hydrological 
regions to evaluate the vulnerability of hydro- 
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logic regions to future climate change. The 
climate change scenarios were obtained from 
two Global Circulation Models (GCMs) and a 
hydrological model was applied to determine 
annual surface runoff, water volume and stor- 
age under future climate conditions (year 2050). 
They showed that the expected climate change 
may have dramatic impacts on the pattern 
and magnitude of runoff, on soil moisture and 
evaporation, and on the aridity level of some of 
Mexico's hydrological zones. Maderey, Jiménez 
and Carrillo (2013) estimated the amount of 
available water over three large basins in Mex- 
ico (i.e., Lerma-Chapala Basin, the Balsas River 
Basin and the Panuco River Basin). Climate 
data were obtained from two GCMs and one 
climate thermodynamic model, and the change 
in water availability was estimated with a water 
balance model for the 2025-2050 period. Their 
results showed that the three climate models 
predict a decrease in available water volumes 
for all basins. 

Rivas-Acosta, Giitrón-De-Los-Reyes and 
Ballinas-González (2010) assessed the climate 
change impact on runoff for three Mexican 
catchments with contrasted climatic conditions: 
The Conchos River Basin (North Mexico), the 
Lerma-Chapala River Basin (Central Mexico) 
and the Grijalva River Basin (South Mexico). 
Rivas-Acosta et al. (2010) considered outputs 
from 23 weighted GCMs under two IPCC 
emission scenarios (A1B and A2). The runoff 
was estimated with a water balance model and 
a vulnerability index was computed. Results 
show that future (2030, 2050, 2100) mean an- 
nual runoff is expected to decrease over the 
three catchments. Tapia, Minjarez and Espinoza 
(2014) evaluated climate change's impacts on 
the water balance of the Yaqui River Basin by 
using climate data from one GCM under two 
IPCC scenarios (A1B and A2) for the 2010-2099 
period. Various runoff behaviors were obtained, 
indicating the possibility of frequent droughts, 
alternating with years of substantially high 
runoff. 

The above studies are based on a low spatial 
resolution of climate models (i.e. GCMs) and 
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the assessments of climate change's impact on 
the catchment's hydrology are estimated using 
a water balance model only. Rainfall-runoff 
models have also been used in the evaluation 
of the climate change impact on water resources 
in Mexican catchments. For instance, Robles- 
Morua, Che, Mayer, and Vivoni (2015) used 
a semi-distributed hydrological model and a 
reservoir optimization algorithm to evaluate 
the hydrological impact of climate change in 
the semi-arid Sonora River Basin (North-West 
Mexico) under the A2 IPPC emission scenario. 
Results showed significantly higher precipita- 
tion in future period (2031-2040) relative to the 
historical period (1990-2000) that would cause 
an increase in reservoir inflows. 

In the present study, the projected change 
in streamflow based on scenarios derived from 
two rainfall-runoff hydrological models is 
explored for a Central Mexican catchment. This 
modeling approach, in which a hydrological 
catchment model is fed with downscaled climate 
scenarios from GCMs, is particularly useful 
for a wide range of impact studies, including 
water resource planning, development and 
management, flood prediction, droughts, water 
quality and hydro-ecology. 

Hydrological models have different degrees 
of complexity and conceptualisation of physi- 
cal processes. They can be classified, consider- 
ing their spatial distribution, as lumped or 
distributed models. In a lumped model, the 
catchment is regarded as a unit. The variables 
and parameters thus represent average values 
for the entire catchment. On the other hand, a 
distributed model takes spatial variations into 
account (e.2., topography, vegetation and soils) 
in all variables and parameters. As previously 
mentioned, hydrological models also differ in 
their representation of physical processes: a 
physically-based model describes the natural 
system using mathematical formulations of 
physical processes while a conceptual model 
is constructed based on physical processes, 
in which physically-based equations are used 
along with semi-empirical equations (Refsgaard, 


1996). Recent research in hydrologic modelling 
tried to take a more physically-based approach 
to understand hydrologic systems” behaviour 
to make better future streamflow predictions 
and to face major challenges in water resource 
management. 

This study aims to evaluate the impact of 
climate change on the hydrology of the Tam- 
paón River Basin in Mexico. In the first step, 
two hydrological models are implemented 
over the basin, and their performances in 
simulating streamflow are compared. The 
selected hydrological models have different 
structural approaches: SWAT is a physically- 
based semi-distributed model (Arnold, Srini- 
vasan, Muttiah, € Williams, 1998) and GR4] 
is a conceptual lumped model (Perrin, 2000). 
Next, climate change impact assessments on the 
catchment's hydrology are evaluated using the 
hydroclimatic model chain illustrated in Figure 
1. This model chain consists of GCM outputs 
dynamically downscaled by a Regional Climate 
Model (RCM) to feed both hydrological models. 
Since RCM outputs usually present consider- 
able biases in climate variables, precluding 
their direct use in hydrological models (Ho, 
Stephenson, Collins, Ferro, $: Brown, 2012), an 
adjustment of the RCM-simulated precipitation 
and temperature is made using a bias correction 
procedure; this allows climate projections to be 
meaningfully translated to the hydrological 
scale (Troin, Velázquez, Caya, zx Brissette, 2015). 
The hydrological simulations are evaluated un- 
der current (1971-2000) and future (2041-2070) 
climates by analyzing hydrological indicators 
for mean, high and low flows. 

The manuscript is organized as follows: first, 
the experimental design is presented, including 
the study area and climate model simulations 
used in this study. Second, a description of the 
bias correction procedure and the hydrological 
models is given. Third, relevant results of the 
catchment-scale hydrological modeling as well 
as an evaluation of the impact of hydrological 
climate change in the Tampaón River Basin are 
analyzed and discussed. Finally, concluding 
remarks close the manuscript. 
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Figure 1. Illustration of the hydroclimatic model chain used 
in this study. 


Experimental design 
Study Catchment 


The Tampaón River Basin is located in the 
East-Central region of Mexico, lying mostly in 
the State of San Luis Potosí (23 373 km?; IMTA, 
2014). It also covers the northern area of the 
States of Guanajuato and Querétaro, and a re- 
gion of the south of the State of Tamaulipas. The 
Tampaón River Basin is a sub-catchment of the 
Panuco River Basin (hydrologic zone 26; Sedue, 
1986) which flows into the Gulf of Mexico. Its 
topographic relief has elevations ranging from 
18 m to 3 500 m with maximal elevation on the 
mountainous western borders (Figure 2). 

The farthest headwater of the Tampaón River 
Basin is the Santa Maria River, which originates 
in the state of Guanajuato, flowing from west 
to east through the Sierra Madre Oriental chain 
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Figure 2. Location of meteorological and discharge gauging 
stations over the Tampaón River Basin. 


mountain. The junction of the Santa Maria River 
with the Verde River makes the start of the Tam- 
paón River. From here, the Tampaón River flows 
to the northeast to join the Gallinas River and 
the Valles River. The Valles River provides the 
largest volume of streamflow to the Tampaón 
River. There are three climatic regions in the 
Tampaón River Basin as a result of topographic 
variation in the Sierra Madre Oriental mountain 
chain (Sedue, 1986): the climate is warm and 
humid in the east, semi-warm and semi-humid 
in the central zone, and dry and temperate in 
the west. 


Observational Datasets 


The daily time series of precipitation (P) and 
minimum and maximum temperatures (Tmin 
and Tmax) were taken from the climatologic 
data base CLICOM (2014), built by the SMN 
(Servicio Meteorológico Nacional de México). 
These data come from seven stations over the 
1971-2000 period (Figure 2, Table 1). The rainy 
season extends from May to October and the dry 
season from November to April (Figure 3). The 
catchment areas mean annual rainfall is about 
1080 mm. The annual average daily temperature 
ranges from 13 to 30 *C with an annual mean 
of 21.5 *C (Figure 4). The coldest and warmest 
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Table 1. Meteorological stations used in this study. 


Station Name 


1 San Luis Potosí 


Villa de Reyes 


Rioverde 


El Salto 


Ballesmi 


2 
3 
Ze Cárdenas 
5 
6 
7 


Xilitla 


600 + r r Y r r r r Y r * ” 


Station 


Mean monthly precipitation (mm month”) 
s 


0 - A 
Jan. Feb. Mar. Apr. May June July Aug. Sep. Oct. Nov. Dec. 


Figure 3. Mean monthly observed precipitation for the 
Tampaón River Basin over the 1971-2000 period. 


months are January and May, respectively. The 
three previously described climatic regions of 
the Tampaón River Basin are well depicted by 
the group of observational data. 

The discharge data comes from the gaug- 
ing station El Pujal (located downstream of 
the Valles River) for the 1971-2000 period; this 
data was obtained from the National Database 
of Surface Water (i.e., Bandas, 2014). Figure 5 
shows the observed mean monthly discharges 
at the gauging station. Two peak flows occur in 
July and September. The decreased discharge in 
August is due to the midsummer drought. Low 
flow values occur from January to May. 


The Climate Model Simulations 


The GCM ensemble simulations used in this 
work are the five members of the Canadian 
Global Climate Model (CGCM3, Scinocca, Me- 
farlane, Lazare, Li, 2 Prammer, 2008) under the 
IPCC SRES-A2 greenhouse gas emission sce- 
nario. Each member of the ensemble had been 
estimated by repeating a climate change experi- 
ment using the GCM several times when only 
the initial conditions were changed by small 
perturbations (Braun, Caya, Frigon, €: Slivitzky, 
2012). Although GCMs are the primary source 
of future climate projections, the spatial resolu- 
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Figure 4. Mean monthly observed temperature for the 
Tampón River Basin over the 1971-2000 period. 
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Figure 5. Observed mean monthly discharges for the 
Tampaón River Basin over the 1971-2000 period. 
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tion of those model outputs is still too coarse to 
be used directly in hydrological impact studies. 
Therefore, the CGCM3 model was dynamically 
downscaled by the Canadian Regional Climate 
Model, version 4.2.3 (CRCM; De Elía € Cóté, 
2010). Daily precipitation, maximum and mini- 
mum temperature series were derived from the 
CRCM simulations for the reference (1971-2000) 
and future (2041-2070) periods. 

The atmospheric component of CGCM3 
is characterized by a spectral model with T47 
horizontal resolution (3.75% x 3.75%) and 31 verti- 
cal layers (Scinocca et al., 2008). The Canadian 
Regional Climate Model derived as an evolution 
of its previous versions (Caya éz Laprise, 1999; 
Laprise, Caya, Erigon, € Paquin, 2003; Plummer 
et al., 2006), which covers the large North Ameri- 
can domain (AMNO); 200 x 192 grid points) with 
a horizontal grid-point spacing of 45-km (true 
at 60 N). The CRCM has largely been used to 
evaluate the impact of climate change on North 
American's catchments (e.2., Muerth et al., 2013; 
Velázquez et al., 2013; Troin et al., 2015), and this 
is the first time that CRCM has been used for an 
impact study on a Mexican catchment. 


Methods 
Bias Correction Procedure 


Climate model simulations are affected by 
biases (í.e., differences between the climate 
model simulations and observations) that, if 
not corrected, can lead to somewhat unrealistic 
streamflow reproduction. Therefore, it is neces- 
sary to perform bias correction on the outputs 
of climate model simulations before using them 
in hydrological models for impact assessment 
(Teutschbein, Wetterhall, €: Seibert, 2011). 
The bias correction method used to adjust 
and downscale CRCM-simulated precipita- 
tion and temperature to the station scale in the 
catchment is the Daily Translation (DT) method 
(Mpelasoka €: Chew, 2009). In DT, different 
correction factors (differences in percentiles 
between observed and climate model simulated 
data during the reference period) are applied 
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to the frequency distribution of projected cli- 
mate model data for the future period. The DT 
method is applied on a monthly basis, and fifty 
percentiles are calculated for each month. The 
corrected temperature and precipitation in the 
reference (ref) period are computed using the 
following equations: 


na a) 0 


re) 
d sim) obs(m,q) sim(m.,g) 


ref ref Led m 1) 
P(corr) "=P; a (2) 
sim(m,q) 


where T (corr) and P (corr) are the bias corrected 
variables, and the subscripts correspond to per- 
centile (9), daily (4) and monthly (1) time steps, 
raw climate simulations (sim) and observations 
(obs). For the future period (fut), the corrected 
precipitation and temperature are obtained 
using; 


T (corr)"" = Tanta) + ELA a Tina) (3) 


Plcorr)'" = pr 


d 


pr 
obs(m,q) 
a) (4) 


sim(m,g) 
Description of the Hydrological Models 


Two hydrological models were selected for 
this study: SWAT and GR4J. The two models 
differ in terms of parameter number, structure 
and physical meaning in their simulations of 
rainfall-runoff processes. Although this study 
provides an opportunity for model comparison 
in a Mexican catchment, both models are also 
applied to the Tampaón River Basin to assess 
the hydrological model's uncertainty in climate 
change impact studies. 

The Soil Water Assessment Tool (SWAT) 
model was developed at the United States 
Department of Agriculture (USDA) by Arnold 
et al. (1998). SWAT is a physically-based semi- 
distributed model that operates the daily time 
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step (Neitsch, Arnold, Kiniry, Williams € King, 
2002). SWAT takes into account the spatial vari- 
ability of the topography, land use, and soil type 
in order to represent the catchment in multiple 
Hydrologic Response Units (HRUs). The input 
variables required to run SWAT are daily pre- 
cipitation and daily maximum and minimum 
air temperatures. The watershed hydrology in 
SWAT is simulated in two steps. The first step 
is the land phase of the hydrologic cycle, which 
calculates the water balance of each HRU to 
provide the amount of water available for each 
sub-basin's main channel at a given time step. 
The second step is channel routing, which 
determines water's progress through the river 
network towards the basin outlet (Neitsch et al., 
2002). A detailed description of model compo- 
nents is presented in Neitsch, Arnold, Kiniry, 
and Williams (2005). 

The GR4J model (which stands for modele 
du Génie Rural 4 4 paramétres Journalier) is a 
daily lumped four-parameter rainfall-runoff 
model developed by Perrin (2000). The input 
variables are daily precipitation and potential 
evapotranspiration (PE). For this study, the PE 
was computed using the formulation (exploit- 
ing extraterrestrial radiation and mean daily 
temperature) proposed by Oudin et al. (2005), 
which has been shown to be as effective for 
rainfall-runoff modelling objectives as more 
complex evapotranspiration formulations. In 
GR4J the hydrology is simulated as follows: 
First, the model subtracts PE from precipitation 
to calculate net rainfall. Then, through intercep- 
tion, a portion of the precipitation goes into the 
production store, where the actual evaporation 
is calculated and percolation occurs. Another 
portion of precipitation goes directly to flow 
routing. The routing part of the structure con- 
sists of two flow components routed by two unit 
hydrographs and a non-linear store. A detailed 
description of the model structure is presented 
in Perrin, Michel and Andreassian (2003). 

The hydrological models” performance is 
evaluated using the Nash-Sutcliffe (NS) coef- 
ficient (Nash éz Sutcliffe, 1970): 


NS= Y (Quin; sl (5) 


Y. (Qui On) 


.and Q. are the observed and 


obs,i sim,i 


where Q 
simulated streamflows at time step í, and N is 
the total number of observations. An efficiency 
of 1 (NS = 1) corresponds to a perfect match of 
modeled discharge to the observed data. 

The second performance criterion is the 
Mean Square Error (MSE) which is computed as: 


N 


(Qu. _ Dore y 
N 


MSE = Lis (6) 


An MSE value of zero corresponds to a 
perfect match between modeled discharge and 
observed data. The MSE is computed by squar- 
ing forecast errors, so itis very sensitive to large 
errors and outliers (e.g., Wilks, 2006). 


Results and Discussion 
Performance of the Hydrological Models 


Table 2 shows the Nash-Sutcliffe (NS) coefficient 
for the calibration (1971-1985) and the validation 
(1986-2000) periods for the two hydrological 
models. Both models present high NS-values, 
ranging from 0.87 to 0.91 in the calibration peri- 
od and from 0.75 to 0.85 in the validation period. 
SWAT performs slightly better than GR4J over 
the calibration and validation periods. 

Figure 6 shows the scatterplot of observed 
versus simulated discharges for SWAT and GR4] 
over the calibration and validation periods. In 
the scatterplot, the more the two data sets agree, 
the more the scatters tend to concentrate in the 
vicinity of the 1:1 line. From Figure 6, it can be 
seen that GR4] better simulates the daily low 
flows, while the daily high flows are better 
predicted by SWAT. As this performance crite- 
rion penalizes errors in the highest streamflow, 
SWAT presents the lowest values of MSE (0.10 
and 0.14) while GR4] presents the highest values 
of MSE (0.14 and 0.23, see Table 2). Figure 6 also 
shows that SWAT tends to slightly overestimate 
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Table 2. Nash-Sutcliffe coefficient (NS) and Mean Square Error (MSE) for the calibration (1971-1985) and validation 
(1986-2000) periods. 


NS MSE 
Hydrological models 
Calibration Validation Calibration Validation 
SWAT 0.91 0.85 0.10 0.14 
GR4] 0.87 0.75 0.14 0.23 
a) SWAT b) GR4] 
10* 104 


hh 
[e] 
a 


hh 
[a] 
mn 


Simulated discharge (més”) 


101 10? 103 10* 
Observed discharge (més) 


Simulated discharge (m%s”!) 


101 10? 103 10* 
Observed discharge (més”) 


Figure 6. Daily discharges as simulated by a) SWAT and b) GR4J for the 1971-2000 period. 


daily streamflow while GR4J slightly underesti- 
mates daily streamflow. 

Figure 7 shows the mean monthly dis- 
charges as simulated by both hydrological 
models driven by observed meteorological data 
over the 1971-2000 period. The most accurate 
simulation of the mean monthly discharges is 
obtained using SWAT. GR4] overestimates the 
mean monthly discharges from November to 
May and underestimates the mean monthly 
discharges from June to September. However, 
both models successfully capture peak flows in 
July and September, and during the midsummer 
drought. 

The performance analysis shows that models 
with different conceptualisation schemes have 
different strengths in simulating the catchment's 
streamflow. The simulation of low flows is 
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Figure 7. Mean monthly discharges for the Tampaón River 
Basin simulated by both models (SWAT and GR4J) using 
observed meteorological data for the 1971-2000 period. 
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challenging_ since water exchange occurs 
through the riverbed and the river may be fed 
by groundwater or may leak to feed the aquifer 
during the low flow period (Pushpalatha, Perrin, 
Le Moine, Mathevet, € Andréassian, 2011). 
On the other hand, the accurate simulation of 
high flows provides further confidence in the 
estimation of extreme hydrologic events, which 
is critical to the evaluation of risk in impact 
studies. 


Future Changes of Precipitation and 
Temperature 


This section presents an analysis of projected 
climate changes over the Tampaón River Basin. 
The predicted changes in precipitation and tem- 
perature in the future climate are determined 
by analysing differences between future (2040- 
2071) and reference (1971-2000) climates derived 
from bias-corrected CRCM simulations. 

Figure 8 shows the mean monthly bias 
corrected climate for the reference period and 
projections for the future period in the three cli- 
matic regions over the catchment. Based on the 


chosen emission scenario, the largest changes 
in mean monthly precipitation are estimated 
during the wet season (June to September). 
For instance, in the dry temperate zone of the 
basin (e.2., climatic region 1), mean monthly 
precipitation during the wet season will vary 
from 55 mm month” to 37 mm month” (Figure 
8a). For the same period, in the semi-humid, 
semi-warm zone of the basin (e.g., climatic re- 
gion 2), precipitation is likely to decrease from 
94 mm month” to 63 mm month” (Figure 8b). 
The most important change in mean monthly 
precipitation is estimated in the warm, humid 
zone of the basin (e.g., climatic region 3) with 
a decrease in mean monthly precipitation from 
323 mm month” to 215 mm month” during the 
wet season (Figure 8c). The projected changes 
in mean monthly temperature are estimated to 
be very similar for the three climatic regions, 
with an increase in annual mean temperature 
of 2.8 *C (Figure 8d to 8f). However, the change 
between the reference and future periods is 
larger for maximum temperature (3.6 *C; not 
shown) than for minimum temperature (2 *C; 
not shown) in all climatic regions. 
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d) Climatic region 1 temp. 
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Figure 8. Mean monthly precipitation and temperature in the reference and future periods from the bias corrected CRCM 


simulations for three climatic regions. 
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Projected climate change signals from the 
bias-corrected CRCM simulations over the 
entire Tampaón River Basin are summarized in 
Figure 9. Future mean monthly temperatures 
are estimated to differ significantly from the 
present time values, with an increase in mean 
temperature between 1.8 and 3.3 "C. The lowest 
temperature rise is found in winter (DJF) and 
the largest temperature increase is estimated to 
occur in summer (JJA). The results also show 
that in the future climate precipitation is very 
likely to increase from October to February and 
to decrease the rest of the year. 


Climate Change Impacts on Streamflow 


Projected climate change's potential impacts on 
hydrology are explored for the Tampaón River 
Basin. Climate-induced streamflow changes are 
inferred by analysing the differences produced 
by SWAT and GR4J when driven by future 
(2041-2070) and reference (1971-2000) climates 
derived from bias-corrected CRCM simulations 
(Figure 10). Based on the considered emission 
scenario, a significant decrease in mean monthly 
streamflow is estimated over the Tampaón 
River Basin as a result of both the increase in 
temperature (which leads to an increase in 
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Figure 9. Mean climate change signals over the Tampaón 
River Basin for the horizon 2070. 
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Figure 10. Mean monthly discharges for the Tampaón River 
Basin as simulated by SWAT and GR4] over the reference 
(1971-2000) and future (2041-2070) periods. 


potential evapotranspiration) and the decrease 
in precipitation. Both models predict reduced 
peak flows in July and September, but the 
models predict different magnitudes. A decrease 
of 309 (286) m? s* is simulated by GR4J in July 
(September), and of 191 (194) m' s* by SWAT for 
the same period. In general, GR4]J predicts larger 
streamflow changes than SWAT all year round. 

The impact of hydrological climate change 
is further evaluated using three hydrologic 
indicators: 


+ Mean monthly streamflow (Qm): the mean 
of all daily values over a given month. 

+ Mean monthly high flow (HF): the mean of 
the maximum streamflow values for a given 
month. 

+ Mean monthly low flow (LF): the mean of 
the minimum streamflow values for a given 
month. 


Climate change's impact on hydrologic 
indicators (Al. ) is expressed as the difference 
between simulated hydrologic indicators from 
the reference (1/%) to the future period (1 ON 
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Figure 11 presents the relative differences 
(AD) in the investigated indicators (Qm, HF and 
LF) over the Tampaón River Basin. GR4J predicts 
a larger decrease than SWAT in all indicators. 
For instance, the relative change in Qm ranges 
from about -30% to -52% during the dry season 
(i.e. December to May) and ranges from about 
-48% to -65% during the wet season (i.e. June 
to September) for SWAT and GR4], respectively 
(Figure 11a). Similarly, the relative change in HF 
ranges from about -29% to -58% during the dry 
season and from about -47% to -70% during the 
wet season for SWAT and GR4J, respectively 
(Figure 11c). Finally, the relative change in LF 
ranges from about -31% to -46% during the dry 
season and from about -48% to -55% during the 
wet season for SWAT and GR4J, respectively 
(Figure 11e). Overall, the median relative dif- 
ferences in the indicators vary from -36% to 
-55% for Qm (Figure 11b), from -34% to -60% 
for HF (Figure 11 d), and from -36% to -49% for 
LF (Figure 11f). These results agree with those 
of Mendoza et al. (1997), which claim that the 
Pánuco River Basin is at risk of becoming a dry 
zone in the future. 


The projected decrease in streamflow over 
the Tampaón River Basin due to changes in 
both temperature and precipitation could have 
a significant impact on water availability for 
economic activities and human consumption. 
Among the economic activities in the catch- 
ment that could be affected by water stress due 
to climate change, are the production of sugar, 
the production of citrus , livestock, industrial 
manufacturing and tourism (Sedue, 1986). In 
addition, the decrease in water availability 
could compromise the water supply of many 
cities in the region. For instance, San Luis Potosí, 
San Miguel de Allende and Celaya will draw 
water from the Santa Maria River (Peña, 2013; 
Conagua, 2012) and it is expected that Monter- 
rey will be supplied by interbasin transfer from 
the Pánuco River Basin (SADM, 2012). 


Conclusions 


An analysis of the projected hydrological chang- 
es is presented for the Tampaón River Basin in 
East-Central Mexico. The selected hydrological 
models are differently conceptualized: SWAT is 
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-20 Es 
3 gs - 
< < -60 EH 
E, 
-80 -80 
Jan. Feb. Mar. Apr. May June July Aug. Sep. Oct. Nov. Dec. SWAT GR4] 
b c) Mean Monthly High Flow (HF) d) HF 
0 
3 201 1 O 
= 404 J Sal E 
< == pa 
-60 + Sl < -60+ 
-80[ 1 1 1 1 1 1 1 1 1 1 1 17 -80 | E 
Jan. Feb. Mar. Apr. May June July Aug. Sep. Oct. Nov. Dec. SWAT GR4] 
e) Mean Monthly Low Flow (LF) f) LF 
20 n 1 1 1 1 1 1 7 1 1 1 1 -20 
S 01 7 840 = 
< 60) 4 <= -60 — E 
-80 1 L L L 1 L L L L 1 L L -80 
Jan. Feb. Mar. Apr. May June July Aug. Sep. Oct. Nov. Dec. SWAT GR4J 


Figure 11. Relative change (%) of the investigated indicators (Qm, HF and LF) for mean monthly flows over the Tampaón River 


Basin. 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VI, núm. 5, septiembre-octubre de 2015, pp. 17-30 


Tecnología y Ciencias del Agua, vol. VI, núm. 5, septiembre-octubre de 2015, pp. 17-30 


Velázquez et al, Hydrological Modeling of the Tampaon River in the Context of Climate Change 


a physically-based semi-distribuited model and 
GR4] is a conceptual lumped model. The study's 
first step was to implement the hydrological 
models over the study catchment. Results show 
that SWAT produces the most accurate simula- 
tion of the observed monthly discharge. GR4J 
overestimates the mean monthly discharge from 
November to May and underestimates the mean 
monthly discharge from June to September. 

The study's second step was to evaluate the 
impact of climate change on the Tampaón River 
Basin's streamflow. Both hydrological models 
estimate an important decrease in mean month- 
ly streamflow as well as in high and low flow 
values for the 2041-2070 period compared with 
the baseline 1971-2000 period. In general, GR4] 
estimates larger streamflow changes than SWAT. 
This result highlights the need to consider the 
uncertainty associated with the hydrological 
models in hydrological climate change impact 
studies. This kind of uncertainty is related to the 
internal computation of hydrological processes, 
such as the estimation of the evapotranspira- 
tion, the soil hydrodynamic formulation and the 
snow model. 

As SWAT has been successfully calibrated 
and validated over that catchment, the imple- 
mentation of the model can be further extended 
to other regions to provide a regional analysis of 
the potential impact of climate change on water 
resources in Mexico. 

The evaluation of climate change impacts 
on streamflow is generally made through a 
model chain that includes GCM outputs for 
a given emission scenario, often dynamically 
downscaled by a RCM and bias corrected with a 
statistical method before using them to a hydro- 
logical model. Each step of the model chain con- 
tributes to the total uncertainty in the estimation 
of future streamflows. The study of Graham, 
Hagemann, Jaun and Beniston (2007) found 
that the most important source of uncertainty 
comes from GCM forcing, which has a larger 
impact on the projected hydrological changes 
than other sources of uncertainty. Blóschl and 
Montanari (2010) pointed out that we can have 
reasonable confidence in predicting hydrologi- 
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cal changes that are driven by air temperature, 
such as snowmelt, than those driven by precipi- 
tation, such as floods. Furthermore, Blóschl and 
Montanari (2010) argued that changing climatic 
conditions have different effects on both the 
weather and hydrology, depending on local fea- 
tures like topography, geomorphology and soils 
of the catchment. Therefore, future work should 
include additional sources of uncertainty, such 
as different GCMs outputs, emission and land 
use scenarios, in order to expand the analysis 
of the uncertainty in the hydrological impact 
studies (e.2., Velázquez € Troin, 2015). 
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